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1 Introduction 

The decay — )■ tt+tt^tt^ is well suited to the study of CP violation and has 
been previously explored by both the BABAR [T] and Belle [2] collaborations. Early 
studies of this mode involved "quasi-two-body" (Q2B) analyses that treated each p 
resonance separately in the decays B° — >■ p^-K^[p^ — t- tt+tt^) and — )• p^7r^(p^ — )■ 
7r^7r°). However, as first pointed out by Snyder and Quinn [3], the use of a full 
time-dependent Dalitz plot (DP) analysis allows sensitivity to the interference effects 
caused by the relative strong and weak phases in the regions where the p"*", p~, 
and p° resonances overlap. This feature makes possible the unambiguous extraction 
of the strong and weak relative phases, and therefore the CP-violating parameter 
a = arg[— V^^VJ^/(\4^\4*^)], where Vqq' are components of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix. A precision measurement of a is of interest 
because it serves to further test the standard model and constrain new physics that 
may contribute to loops in diagrams. 

In this paper, we summarize an extensive reoptimization of an earlier BABAR 
analysis. We use the full "on-resonance" BABAR dataset of approximately 431 fb~ 
collected at the T{AS) resonance (an increase of 25% in the number of B decays) 
and include a number of improvements to both reconstruction and selection. Among 
these are improved charged-particle tracking, improved particle identification (PID), 
and a reoptimized multivariate discriminator (used both for event selection and as a 
variable in the final fit). 

^Speaker on behalf of the BaBAR Collaboration 

^Throughout this paper, whenever a mode is given, the charge conjugate is also implied unless 
indicated otherwise. 



2 Reconstruction and Event Selection 



The data used in this analysis were collected with the BaBar detector at the PEP-II 
asymmetric-energy e~^e~ storage ring at SLAC. Collisions occur at the T{4S) reso- 
nance energy {y/s ~ 10.58 GeV/c^), which frequently decays to BE pairs. We fully 
reconstruct the decay of one B (.63^) and use the decay of the other B (-Btag) to 
determine the flavor of 5tag at the time of its decay. Due to the asymmetric energies 
of the and beams, the e+e^ centcr-of-mass (CM) has a boost of f3'~f ^ 0.56 
in the laboratory frame. The time-dependence of our analysis is measured using the 
distance along the beam axis between the B^^^ and Stag decay vertices to calculate 
the time At between the two decays. 

Pairs of oppositely charged tracks are combined with 7r° — >■ 77 candidates to 
construct B^ — ?■ tt^tt^tt^ candidates. The kinematics of B meson decays that are 
fully reconstructed at BaBar can be characterized by two variables: rriEs and AE. 
The beam-energy-substituted mass uies is the invariant mass of the reconstructed B 
candidate calculated under the assumption that its energy in the e+e~ CM frame is 

half the total beam energy. We define m-Es — y^[(s/2 +pi ■ pB)/Eif' — where 
^/s is the total beam energy in the e~^e^ CM frame, (Ei,pi) is the four-momentum of 
the e+e~ system in the laboratory frame, and pb is the i?-candidate momentum in 
the laboratory frame. The second kinematic variable is defined by AE = Eq — \\/s, 
where E^ is the measured energy of the B candidate in the e+e" CM frame. We 
apply loose selection criteria using these variables and include them as inputs to the 
final fit. 

Basic selection criteria are applied using quantities such as photon lateral mo- 
ments (for the 7r° 77 candidate), energy deposits, and track geometry parameters. 
Additionally, we use PID information to require that the candidates be consistent 
with the pion hypothesis. We apply a loose selection criterion using At and include 
it as a variable in the final fit. 

A further selection criterion is applied using a multivariate neural network (NN) 
discriminator. The discriminator serves to distinguish signal-like events (which have a 
more spherical topology) from qq continuum (g = m, d, s, c) background events which 
have a more coUimated event shape. We train the discriminator using signal B^ — > 
pTT Monte Carlo (MC) and data collected below the B'^B threshold (to represent 
continuum background). A loose selection criterion is applied using the NN and we 
include it variable in our final fit. 

While one would typically parameterize a Dalitz plot using the squared invariant 
masses of two pairs of daughter particles, practical considerations lead us to use a 
square Dalitz plot (SDP) parameterization in which the kinematically allowed region 
of DP phase space is mapped onto a unit square. The square Dalitz plot coordinates 
are m', which depends on tt+tt" invariant mass, and 6*', which depends on the p° 
helicity angle 6q. 
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3 Maximum Likelihood Fit 



We perform an unbinned extended maximum likelihood fit in order to extract event 
yields and physics parameters. The input variables are mss, AE", the NN output, 
and the three time-dependent-SDP variables m', 9' , and At. We also use (Ta* (the 
per-event uncertainty on At) as a scale factor in the signal At resolution function. 
The likelihood function used in the fit consists of separate components for signal, 
continuum background, charged B backgrounds, and neutral B backgrounds. The 
signal component is subdivided into correctly reconstructed and misreconstructed 
components. 

A probability density function (PDF) is associated with the distribution of each 
fit variable in each component of the likelihood function. Fixed and initial parameter 
values for these PDFs are obtained from fits to fully simulated MC (for signal B 
background components) and either data collected below the B^B^ threshold or a 
lower sideband in tties (for continuum). 

We parameterize our signal PDF using 27 real-valued U and / coefficients, defined 
in terms of B^ — t- p+'^'Oyr^'+'O and B — > p+'^'Oyr^'+'O decay amplitudes {A'^~'^ and 



A^' '°, respectively) as = \A''\^ ±\A''\\ U^^^^^"''^ = Re(Im) 



= Im 



A A 



/l^ = Re 



and J,™ = Im 



A" A"* ± A A 
A'" A"* + A" A 



A A''* -A A^* 

where k G (+, — , 0) and G (H — , +0, —0). These coefficients provide an alternative 
parameterization to tree and penguin amplitudes (as well as a) or to the amplitudes 
A'^ and A'^ The U and / parameters can also be directly related to the Q2B 
parameters {C, S, AC, AS, Ap-,,, Cqq, Sqq, and /oo) often used in CP- violation analy- 
ses [5]. 

We include the p(1450) in the final fit with an assumption that the relative mag- 
nitudes and phases between the three p(1450) resonances are the same as for the 
p(770). Whereas there is reasonable motivation for this assumption in the case of the 
p(1450) since the p(770) and p(1450) have the same quantum numbers, the p(1700) 
does not share these quantum numbers (£ = 2 instead of 0). Since the p(1700) is not 
expected to have a large contribution to the decay rate, we excluded the p(1700) from 
the fit and associate a systematic uncertainty with this omission. We also associate 
systematic uncertainties with the expected numbers of B background events (which 
are fixed in the fit), the masses and widths used in our p lineshapes, possible contri- 
butions from uniform backgrounds, and various other small contributions. While the 
most significant systematic uncertainty is that associated with the exclusion of the 
p(1700) from the nominal fit model, even this contribution is found to be small and 
the sensitivity remains dominated by statistical uncertainties. 
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4 Results 



From an on-resonance dataset containing 53,084 candidates, the fit extracts 2,940±100 
signal events and 46,750±220 continuum events. Figure [T] contains overlaid vrvr invari- 
ant mass plots of the data used in the final fit and parameterized MC generated using 
the results of the final fit. A study of the U and / parameters (see Sec. [s]) exhibits 
negligible bias in their extraction, and good robustness in the presence of statistical 
fluctuations. 




m(it* jt ) (GeV/c') m(ii* it°) (GeV/c') m(7t- (GeV/c^) 

Figure 1: Overlay plots of on-resonance data (points with error bars) and parameter- 
ized MC generated from the final fit results (red line) with 10 times the number of 
events in data. The MC histograms are scaled to have the same integral as the data 
histograms onto which they are overlaid. A tight cut is applied to the NN output to 
enhance the signal component. 

The complete final results from the extended maximum likelihood fit, including 
statistical and systematic uncertainties, are provided in the Physical Review D article 
associated with this analysis (currently in preparation). We find the extracted U and 
/ parameter values as well as the extracted Q2B parameter values to be consistent 
with the previous results from BaBar and Belle. The sensitivity of the present analysis 
is much improved over the previous BaBar analysis, with an average ratio of statistical 
uncertainties on U and / parameters relative to the previous analysis of 0.47. This 
is a larger increase in sensitivity than can be explained by the 25% increase in the 
size of the data sample and it may be attributed to the many improvements made 
in this analysis. Similarly, the average ratio of the statistical uncertainty on the 
Q2B parameters in this analysis relative to the previous BaBar analysis is 0.61. The 
final fit results for the Q2B parameters are provided in Table [1} A study of the 
Q2B parameters (see Sec. |5| exhibits negligible bias in their extraction, and good 
robustness in the presence of statistical fiuctuations. 

The parameters Apn and C can be transformed into the direct CP-violation pa- 
rameters Ap~ and as described in Ref. [I]. We extract the central values and 
uncertainties for these parameters using a minimization in the two-dimensional 
plane corresponding to A't^ vs. A'^J'. From this two-dimensional scan, we find 



4 




Figure 2: Combined statistical and systematic two-dimensional likelihood scan of 
vs. with 68.3%, 95.5%, and 99.7% confidence-level contours (Ax^ = 

{2.30, 6.18, 11.83}). The yellow point indicates the central value. 
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0.016 


0.059 
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0.19 
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0.234 
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-0.37 


0.34 


0.20 


S 


0.053 


0.081 


0.034 


/oo 


0.092 


0.011 


0.008 



Table 1: Quasi- two-body parameter values and uncertainties corresponding to the fit 
to the complete on- resonance dataset. 

A+- = Omttfe ± 0.04 and A-+ = -0.12 ± 0.08ll];[Jt A plot of this scan is pro- 
vided in Fig. [2j The origin corresponds to no direct CP violation and lies on the 
96.0% confidence level contour {Ax^ = 6.42). The corresponding p value for the 
hypothesis of no direct CP violation is 4.0%. 

4.1 cx Scan Results 

In order to extract likely values of a from the U and / parameters obtained in our final 
fit, we perform a scan of a from 0° to 180°. At each scan point, a x^-minimization 
fit is performed using the combined statistical and systematic covariance matrices 
from our nominal fit. As the scan proceeds, a minimum value is extracted from 
the fit at each value of a. We convert these values to "Z"' values by calculating 
the probability of each value according to = /(x; l)dx, where a is the 
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difference between the at the current scan point and the minimum for all the 
scan points, and /(x; 1) is a distribution with one degree of freedom. The variable 
"Z"' corresponds to what is commonly referred to as "1— Confidence Level" (1— C.L.). 

Following the methods employed in Belle's 2007 — )• pvr analysis P] and de- 
scribed in [4j, we perform a further a scan that makes use of measurements from 
the charged decays — )■ p^'^vr^'^. Amplitudes for these modes can be related to 
amplitudes in the neutral B modes due to isospin relations. These relations result in 
four constraint equations while introducing only two new free parameters in the fit 
(which arise from the unknown relative magnitude and phase of the charged- -B and 
neutral- -B decay amplitudes). 

Graphs of the x^ values from our final a scans with isospin constraints (solid red) 
and without isospin constraints (dashed black) are provided in the left plot of Fig. |3} 
The corresponding S distributions are given in the right plot of Fig. [3j Importantly, 
our robustness studies (see Sec. |5]) indicate that the E scan is not robust with our 
current sample size (or those available to the previous BaBar and Belle analyses) and 
cannot be interpreted in terms of Gaussian statistics. 




Figure 3: Left: Scans of minimum values as a function of a. Note that the 
origin on the vertical scale is suppressed. Right: Scans of U values as a function of 
a. The upper and lower horizontal dashed lines correspond to = 0.05 and 0.32, 
respectively. All scans are based on the fit to the full on-resonance dataset and include 
contributions from both statistical and systematic uncertainties. 



5 Robustness Studies 

An important component of this analysis is a set of studies which assess the robustness 
with which the fit framework extracts statistically accurate values and uncertainties 
for the U and / parameters, the Q2B parameters, and a by employing 25 MC simu- 
lated samples generated with a parameterized detector simulation and with signal and 
background contributions corresponding to those expected in the full on-resonance 
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dataset. The samples are simulated using physical U and / parameters generated 
based on specific tree and penguin amplitudes and a = 89° (approximately the world 
average). Each simulated dataset is generated with the same parameter values, but a 
different random-number seed. By examining the results of fits to each of these sim- 
ulated datasets, we assess the robustness of the fits. Comparisons of the extracted U 
and / and Q2B parameter values with the generated values find that all parameters 
are robustly extracted with negligible bias and well estimated uncertainties. More 
significant are the results of the a robustness study. 

A one-dimensional likelihood scan of a is performed using the results of the fits 
to each of the 25 MC samples. For 8 of the 25 scans, the extracted value of a lies 
more than 3a from the generated value. Examining the individual a scans reveals 
three distinct solutions for a that tend to be favored (including the generated value 
of 89°) and each scan tends to include at least one secondary peak in addition to 
the primary peak. The left plot in Figure |4] illustrates the three solutions for a by 
providing the sum of 25 normalized Gaussians with means and widths determined 
by the peak positions and symmetric errors extracted from the 25 a scans. Because 
the errors are not truly Gaussian, the plot provides an incomplete picture of the scan 
results. A better illustration is provided by the right plot in Figure |4| which displays 
the total E distribution obtained by summing all 25 a scans after normalizing each 
to the same area. The total distribution is scaled so that it peaks at 1. The final PDF 
closely resembles that obtained by naively summing Gaussian distributions, though it 
exhibits more fine features. Again, the distribution indicates three distinct solutions 
for a, with the generated value of 89° being favored. At the la level (17=0.32), the 
total scan distribution allows both the central and left peak. The presence of these 
secondary solutions indicates that with the current signal sample size and background 
levels, there is still a significant possibility that the favored value of a in a particular 
scan will correspond to a secondary solution. 

6 Conclusions 

We have performed a time-dependent Dalitz plot analysis of the mode 5° — (p7r)° in 
which we extract 26 U and / parameter values describing the physics involved, as well 
as their full statistical and systematic covariance matrices in an extended unbinned 
maximum likelihood fit to the full BaBar dataset. From these fit results, we extract 
standard quasi-two-body parameters with values given in Table [Tj These Q2B values 
are consistent with the results of the 2007 BaBar [T] and Belle analyses, but exhibit 
significantly increased sensitivity. We also perform a two-dimensional likelihood scan 
of the direct CP-violation asymmetry parameters for — )■ p^vr^ decays, finding the 
change in between the minimum and the origin (corresponding to no direct CP- 
violation) to be Ax^ = 6.42 (approximately 2a). Finally, we perform one-dimensional 
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Figure 4: Left: Plot of the sum of 25 normalized Gaussians based on peaks from 
scans of parameterized MC generated with a = 89°. Also plotted are the individual 
Gaussians which contribute to the total PDF. Right: Plot of the total 2J distribu- 
tion obtained by summing up all 25 E scans of parameterized MC generated after 
normalizing their areas to 1. Also plotted are the individual E scans. 

likelihood-scans of the unitarity angle a both with and without isospin constraints 
from other modes (see Fig. [3]). 

Notably, we also perform a series of robustness studies in order to determine 
how reliably our fit framework extracts the actual value of physics parameters. The 
studies reveal that with our current signal sample size and background suppression, 
we can reliably extract the U and I parameters as well as the Q2B parameters, 
but the extraction of a in 5° — t- pvr is not statistically robust. This result has 
consequences not only for this analysis, but earlier BaBar and Belle published results 
as well. Namely, it calls into question the reliability of the a values quoted in previous 
analyses and indicates that with the current sensitivity, we should not treat the a 
scan as a straightforward measurement of a. This analysis would benefit greatly from 
increased sample sizes available at high-luminosity experiments such as Belle II. 
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